(N 



£ 



^o 



o 



The High Density Effects in Heavy Quark 
Production at pA Colliders 



"j- AT Aomin Pilhn M „„,» V P CZn^rnloyPQ l 2'** 

1 Instituto de Fisica e Matemdtica, Univ. Federal de Pelotas 
Caixa Postal 354, 96010-090, Pelotas, RS, BRAZIL 



2 Universidade Estadual do Rio Grande do Sul - UERGS, 
Borges de Medeiros, 1501, 90119-900, Porto Alegre, RS, BRAZIL 



pq . Abstract: In this paper we investigate the role of the high density 

effects in the heavy quark production cross section in pA processes at RHIC 

^ and LHC We use, as initial condition, a gluon distribution consistent with 

fixed target nuclear data and the Glauber-Mueller approach to describe the 



high density effects. We show that this process can be used as a probe of 
the presence of the high density effects. Moreover, we include these effects in 
the calculation of the heavy quark production in AA collisions, verifying that 
they cannot be disregarded both in the estimates of quarkonium suppression 

and in the initial conditions of the quark-gluon plasma. 
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The relativistic collider facilities RHIC and LHC will for the first time 
provide the opportunity to systematically study the physics of hot and ultra- 
dense matter in hadron-nucleus (pA) and nucleus-nucleus (AB) collisions at 
energies that are orders of magnitude larger than the energies of the current 
accelerators. The systematic study of pA collisions at the same energies is 
essential to gain insight into the structure of the dense medium effects. Such 
effects, as the energy loss and shadowing, are absent or small in pp collisions, 
but become increasingly prominent in pA collisions, and are of major impor- 
tance in AA reactions. By comparing pA and AA reactions involving very 
heavy nuclei, one may be able to distinguish basic hadronic effects that dom- 
inate the dynamics in pA collisions, from a quark-gluon formation predicted 
to occur in heavy ion AA collisions. To gain insight into the underlying 
hadronic processes, one has to study collisions that are expected to not lead 
to a QGP formation. Once the physics of "QCD at high densities" is better 
understood, the mechanisms of quark-gluon plasma formation and related 
collective phenomena in heavy ion collisions could be disentangled from the 
basic hadronic effects. 

One of the nuclear medium effects is the nuclear shadowing, which is 
the modification of the target parton distributions so that xq A (x,Q 2 ) < 
Axq N (x,Q 2 ), as expected from a superposition of pp interactions. The cur- 
rent experimental data presenting nuclear shadowing can be described rea- 
sonably using the DGLAP evolution equations |L| with adjusted initial par- 
ton distributions ||. However, this parameterization does not include the 
dynamical saturation effects of the parton distributions predicted at small x 
and large A ||. In this kinematical regime, the density of quarks and gluons 
becomes very high and the processes of interaction and recombination be- 
tween partons, not present in the DGLAP evolution, should be considered. 
Recently, we have proposed a procedure to improve the nuclear parton distri- 
butions and include the perturbative high density effects in the predictions of 
the inclusive observables in eA processes [Q. There we have indicate that the 
high density effects becomes very important mainly for the LHC kinematic 
region and that these effects could not be disregarded in the calculations of 
the observables and signatures of the QGP. Here we analyze in detail the 
heavy quark production in pA processes in order to investigate the presence 
and magnitude of the high density effects (For a related discussion of heavy 
quark production see Refs. || ||). Our analyzes is motivated by the fact 
that the heavy quark production at RHIC and LHC energies is dominated by 

1 



initial state gluons. Thus, the pA process becomes a valuable search of infor- 
mation about the gluon distribution in the nucleus. Since the probability for 
making a heavy quark pair is proportional to square of gluon distribution, 
any depletion in number of gluons will make a significant difference in the 
number of heavy quark pairs produced. 

Lets start from a brief review of the high density scenario. At high den- 
sities we expect the limitation on the maximum phase-space parton density 
that can be reached in the hadron wavefunction (parton saturation) and very 
high values of the QCD field strength squared F 2 oc l/a s [|J, possibly form- 
ing a Color Glass Condensate || , which is characterized by a bulk momentum 
scale Q s . If this saturation scale is larger than the QCD scale Aqcd, then 
this system can be studied using weak coupling methods. At present, there 
are different approaches to describe hdQCD effects, all of them including 
nonlinear terms in parton evolution equations ||. In fact, these approaches 
sum powers of the probability of gluon-gluon interaction inside of the par- 
ton cascade. As a summary, we have that the different approaches match 
the DLA limit of the DGLAP parton evolution in the low parton density 
limit and match the Glauber-Mueller approach in the transition regime of 
low density to the high density limit. At the energies of interest, we believe 
that the predictions of the Glauber-Mueller approach can be taken as good 
approximation to the description of high density effects. As this approach is 
a common limit to the distinct approaches, we intend to obtain predictions 
which are not model dependent. 

The main properties of the nuclear gluon distribution xG A (x, Q 2 ) in the 
Glauber-Mueller approach were extensively discussed in Ref. [[UJ. For com- 
pleteness, we now present a qualitative discussion of the main properties of 
this approach. In the nucleus rest frame we can consider the interaction 
between a virtual colorless hard probe and the nucleus via the gluon pair 
(gg) component of the virtual probe. The interaction of the dipole with the 
color field of the nucleus will clearly depend on the size of the dipole. If the 
separation of the gg pair is very small (smaller than the mean separation of 
the partons), the color field of the dipole will be effectively screened and the 
nucleus will be essentially 'transparent' to the dipole. At large dipole sizes, 
the color field of the dipole is large and it interacts strongly with the target 
and is sensitive both to its structure and size. More generally, when the 
parton density is such that the nucleus becomes 'black' and the interaction 
probability is unity, the dipole cross section saturates and the gluon distribu- 



tion becomes proportional to the virtuality of the probe Q 2 . In the infinite 
momentum frame, this picture is equivalent to a situation in which the in- 
dividual partons become so close that they have a significant probability of 
interacting with each other before interaction with the probe. Such interac- 
tions lead, for instance, to two — > one branchings and hence a reduction in 



the gluon distribution. These properties were considered in Ref. ||10| , where 
the rescatterings of the gluon pair inside the nucleus were estimated using 
the Glauber-Mueller approach, resulting that the nuclear gluon distribution 
is given by flO 
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xGa (x,Q 2 ) = ^ ~— ?—±{C + ln(K G (2<,ri)) + E 1 (K G (x',ri))} (1) 



where C is the Euler constant, E\ is the exponential function, the function 
K G (x,r 2 ) = (3a s A/2R\) nr 2 xGn(x, -^), A is the number of nucleons in a 

nucleus and R\ is the mean nuclear radius. One of the shortcomings of this 
approach is that the Glauber-Mueller formula disregards any nuclear effect in 
the nonperturbative initial condition for the gluon distribution. Furthermore, 
the medium effects present at larger values of x, the antishadowing and the 
EMC effect, are also disregarded in this approach. In order to improve the 
Glauber-Mueller approach, it was recently proposed in Ref. f| a modification 
in the expression ([TJ) that includes the full DGLAP kernel in parton evolution. 
Basically, we propose to calculate the nuclear gluon distribution using the 
following procedure: 

xG A {x, Q 2 ) = (l/A)xG A (x, Q 2 ) [GM] - (l/A)xG A (x, Q 2 ) [DLA] 

+ xG A (x,Q 2 )[EKS] , (2) 

where xG A (x, Q 2 )[GM] represents the Glauber-Mueller nuclear gluon dis- 
tribution given by the expression (|l|) and xG A (x, Q 2 )[DLA] is the DGLAP 
(DLA) prediction for the nuclear gluon distribution, which correspond to 
the first term of expression ([[]) when expanded in powers of Kq. The last 
term in expression (g), xG A (x,Q 2 )[EKS], is the gluon distribution solu- 
tion of the DGLAP equation as proposed by Eskola, Kolhinen and Sal- 
gado(EKS) ||, where the initial conditions of the parton evolution are cho- 
sen in such a way to describe the nuclear effects in DIS and Drell-Yan fixed 
nuclear target data. As it was discussed in detail in Ref.Q, the parame- 
terization of the nuclear effects in the gluon distribution, given by the ratio 



Rq(x, Q 2 ) = xGa{x, Q 2 )/xGn(x, Q 2 ) in the EKS procedure, does not include 
the perturbative high density effects at small values of x. Thus, the proce- 
dure presented in equation (0) includes the full DGLAP evolution equation 
in all kinematic region, take into account the nuclear affect in the present 
fixed target data and includes the high density effects in the parton evolution 
at small x in the perturbative regime. 

To illustrate our results, in Fig. |T] we compare our predictions for the 
x - dependence of the ratio Rg = xGa/xGn for Q 2 = 15 GeV 2 with the 
EKS predictions. The results are shown for two typical nuclei of interest 
in nuclear collisions. Also the corresponding kinematic region which will 
be explored in the RHIC and LHC colliders are presented. Two important 
features of the high density effects can be seen in Fig. [I]: the magnitude 
of the effects and the saturation behavior. The suppression of the gluon 
distribution in respect to the EKS set is 6% for A=40 (Ca) and 11% for 
A=208 (Pb) in the lower limit of the RHIC kinematic range. When the lower 
limit of the LHC kinematic range is concerned, the reduction of the gluon 
distribution is 31% for A=40 and 50% for A=208. This strong effect also 
modify the saturation of the ratio predicted in the EKS parameterization at 
low x. Basically, in the EKS approach the saturation of the ratio is assumed 
in the nonperturbative initial condition as a constraint for the behavior of the 
nuclear parton distributions in the small x limit. This behavior is preserved 
by the DGLAP evolution, since these equations reduce to the DLA limit at 
low x in the nucleon and nuclear case, keeping the ratio constant. When the 
high density effects are considered in the dynamics, a large modification of 
the gluon distribution is predicted, which is amplified in nuclear processes 
since the nuclear medium amplifies the effects associated to the high parton 
density. Therefore, the high density approaches predict larger effects in the 
nuclear gluon distribution than the nucleon one, independent of the initial 
parton distributions used. This expectation is present in the behavior of 
Rq(AG), which turn to be much smaller then R G (EKS) in low x region. 
As the heavy quark production is dominated by the fusion of gluons, we 
expect that these effects will strongly modify the cc and bb production in pA 
processes in the high energy limit. 

Now let us consider the heavy quark production. Perturbative QCD cal- 
culations of heavy quark production at leading order have long been available. 
For high energies and at leading order (LO), the production is dominated by 
gluon fusion, with charm and bottom quarks produced basically in the pro- 



cess gg — > QQ with Q = c and b. The LO cross section for a proton p 
colliding with a nucleus A is then: 



0>a = dr dxF , xG p (xi, [/) xGa(%2, fJ 2 ) x 

%-,QQ(^i^2,m Q ,^ 2 ) , (3) 



where xp = x\ — %2, t = M 2 / s and M 2 is the invariant mass of the virtual 
gluon in the subprocess. xGi the gluon distributions evaluated at momentum 
fraction x and momentum scale fi 2 . For practical purpose, we take fi 2 = 
Am 2 for charm production and fi 2 = ml for bottom production. In order 
to investigate the medium dependence of the heavy quark production cross 
section, we will follow the usual procedure used to describe the experimental 
data on nuclear effects in the hadronic quarkonium production ||12| , where 
the atomic mass number A dependence is parameterized by <j v a = o" p at x A a . 
Here a p A and a P N are the particle production cross sections in proton-nucleus 
and proton-nucleon interactions, respectively. If the particle production is 
not modified by the presence of nuclear matter, then a = 1. A number 
of experiments have measured a less than linear A dependence for various 
processes of production, which indicates that the medium effects cannot be 
disregarded. As our focus is to analyze the influence of the high density 
effects in the heavy quark production when compared to the usual DGLAP- 
EKS predictions, we do not consider the contributions of energy loss and 
the intrinsic heavy-quark components for the non-linear A dependence of the 
cross sections (For a discussion of these effects in charmonium production 
see Ref. |H|). In this paper we will evaluate the heavy quark production at 
leading order and will use the quark masses m c = 1.2 GeV and m& = 4.75 
GeV. To estimate the modification of heavy quark production cross section 
due to the high density effects, we calculate the effective exponent a, which 
is given by 

a = (ln^-1 'in a) . (4) 

We calculate the above expression considering, as an input in the nuclear 
cross sections, the parameterization of the nuclear effects in the gluon dis- 
tributions, which is given by xGa(xi, fJ> 2 ) = Rc{. x i^ 2 ) %G n(xi, A* 2 ) , with Rq 



presented in Fig. [J] and xGn taken from the GRV94(LO) set [[HJ]. Two 
comments are in order here. First, since we are interested in the behavior 
of the effective exponent, we may expect that our results will not be mod- 
ified by the NLO corrections (the K- factor) fLlf . Second, we are assuming 
that the collinear factorization is still valid in the kinematic region consid- 
ered. This is a strong assumption which should be tested by the comparison 
of our results with a calculation of this process considering the color glass 



condensate formalism, similarly to made in Ref. [18| for photoproduction. 
The agreement between the results will allow us to use our approach as a 
simplified method of inclusion of the high density effects in the nuclear cross 
sections. Furthermore, the comparison between the results of our approach 
and the color glass predictions for the gluon minijet production [Tj| will be 



other important cross check of our approach. Work in this direction is in 
progress. 

In order to gain insight into the amount of the high density effects in heavy 
quark production at future colliders, we will consider the process calculated 
at the proposed energies for RHIC and LHC These energies are a/s = 200 
GeV per nucleon pair for Au+Au collisions and yfs = 350 GeV for pA process 
at RHIC. For LHC, \/s = 5500 GeV per nucleon pair for Pb+Pb collisions 
and \/s = 8800 GeV for pA processes |20). In Fig. |] we present the effective 



exponent a as a function of the cm. energy s 1 ^ 2 for two different nucleus. 
We verify that the high density effects are sizeable at high values of energy 
even for a small value of the mass number (A = 40). The general behavior 
of the exponent can be understood as follows. When the integrations in Eq. 
(H) are taken, the nuclear gluon distribution are evaluated in the x<i interval 
given by r < ^2 < \fr. Thus, when the energy grows, the x<i interval goes 
to the small x region. For charm production, for example, the antishadowing 
region in Fig. [L] dominates the integration for energies smaller than 80 GeV. 
For bigger energies, the shadowing region dominates and the exponent is 
smaller than 1. In bottom production, the x^ interval is dislocated to larger 
values of x, which implies that the antishadowing region of Fig. [1] dominates 
the production at energies in the interval 100-200 GeV and a is close to 
one even for RHIC energies. For bigger values of energy, the suppression is 
sizeable and the effective exponent is smaller than 1. In general grounds we 
have that the saturation of the exponent observed in the EKS prediction is 
associated to the behavior of the ratio Rg at small x, while the presence of 
the high density effects (nonsaturation of Rq) implies a large reduction of the 

6 



pA cross section when compared with DGLAP-EKS description of nuclear 
effects. Therefore, we believe that the analyzes of the effective exponent in 
pA processes can be useful to evidence the high density effects. 

In Fig. H we present the exponent a for RHIC and LHC energies as a 
function of xp = x\ — x^- The exponent was obtained from the expression 



*M = \ln ^ r^\ lnA\, (5) 




where the differential cross section was obtaining integrating equation (|3]) in 
the invariant mass of the virtual gluons from the heavy quark threshold to the 

2 

energy squared s, which corresponds to -# < r < 1. For values of xf close 
to zero, X2 ~ X\ and the r integration corresponds to a large interval in both 
X2 and x\. For xf close to one, both X2 and x\ interval are much smaller, 
with x\ close to one and X2 close to zero. The lower limit of r corresponds 
also to the lower X2 limit. Considering the kinematics discussed above and 
the properties of the Fig. [L], all the features of Fig. |3] can be understood. For 
the p — Pb processes at RHIC (\/s = 350 GeV) with xp close to zero, the 
exponent is close to one. This occurs because the r integration corresponds to 
a large interval in X2, including contribution from the antishadowing region 
of Fig. |TJ. For charm production, the threshold is low enough to include 
the contribution from the small X2 region and the exponent a is smaller 
than one. For the bottom production, the low X2 contribution is small and 
the exponent is bigger than one. Since high density effects are small in the 
intermediate X2 region, the EKS and AG prediction are close to each other. 
For xp close to one, the low X2 behavior of the gluon distribution dominates 
the process. Since the EKS parameterization predicts that the ratio saturates 
at this region, the exponent tend to a constant, while the high density effects 
predict an exponent that decrease monotonically. These effects are much 
stronger in the charm production due to the lower threshold. Moreover, 
all effects discussed above are amplified for LHC energy (y/s = 8.8 TeV), 
even for xf close to zero. Due to the large energy, the lower limit of the 
r integration is very close to zero. Then, the low X2 contribution is much 
bigger and the high density effects are much more important. For xf close 
to one, the EKS prediction indicates a constant value for a, both for charm 
and bottom production. This is a direct consequence of the saturation of 
the ratio Rq in the low x 2 region. As far as the high density effects are 

7 



considered, the exponent is strongly reduced. As xf tend to one, a goes to 
0.6 for charm production and 0.7 for bottom production. 

The strong modification of the heavy quark production in pA process indi- 
cate that high density effects on the gluon distribution will play an important 
role in the calculation of the initial conditions and signatures of the quark- 
gluon plasma in heavy ion collisions. To illustrate this point, we present in 
Fig. | the exponent (3 calculated by the expression 

0= (ln^± I ' lnA\ , (6) 

for the heavy quark production in Pb-Pb collisions as a function of ■y/s. The 
exponent present a similar feature to the exponent a shown in Fig. |2|. For 
energies close and smaller than 100 GeV, the most important contribution 
comes from the antishadowing region of Fig. [J], and the exponent is bigger 
than two. For larger energies, the small x region gives the more important 
contribution. When compared to the EKS prediction, the high density effects 
strongly reduces the production cross sections. At RHIC energy (-^/i = 
200GeV), the exponent /3 is reduced in 8% for cc and 4% for bb. At LHC 
energy (s/s = 5500), the relative reduction is 42% for cc and 28% for bb. 
As we can see, the high density effects will strongly reduce the production 
cross sections in AA processes and should be taken into account to calculate 
the cross section of the hard QCD process in the first stage of the heavy ion 
collisions. 

In brief, a systematic measurement of gluon shadowing is of fundamental 
interest in understanding of the parton structure of nuclei as well as in the 
field of minijet production that determines the total entropy produced at 
RHIC and higher energies. In this paper we addressed the heavy quark pro- 
duction in pA processes as a search to verify the presence and estimate the 
magnitude of the high density effects and, consequently, fix the behavior of 
the gluon distribution at large energies. In this sense our results can be con- 
sidered as complementary to the predictions made by Wang and Gyulassy 



21j| , which have studied the sensitivity of single-particle inclusive spectra 
in nuclear collisions to gluon shadowing and jet quenching, and have sug- 
gested the investigation of the pA collisions at the same energy in order to 
disentangle both effects. Moreover, the study of heavy quark production is 
important to estimate the suppression of this process associated to the high 



density effects, which is necessary for a reliable calculation of quarkonium 
production and its suppression in a quark-gluon plasma. As a summary of 
our results we have calculated the cross sections for heavy quark production 
in pA processes and estimated the energy and xp dependencies of the effective 
exponents that parameterize the medium effects. Our results demonstrate 
that these effects are large even at small nuclei and that a systematic ex- 
perimental analyzes could discriminate between the high density regime and 
the predictions from the linear regime. Furthermore, we have extended our 
analyzes for heavy quark production in heavy ion collisions and estimated 
the contribution of the high density effects in this process. Our results have 
important implications in the signatures of the QCD phase transition, spe- 
cially in the suppression of quarkonium production since our results indicate 
that the quarkonium production rate should be strongly modified by the 
presence of the high density effects. Consequently, if the analyzes of the pA 
process demonstrate the presence of the these effects in the kinematic re- 
gions of RHIC and LHC, we have that the current estimates of quarkonium 
suppression in AA processes should be completely reanalyzed. 

Acknowledgments 

The authors acknowledge helpful discussions with M. B. Gay Ducati, J. 
Jalilian-Marian, M. V. T. Machado and R. Venugopalan. VPG thanks the 
Brookhaven National Laboratory for its hospitality at a preliminary stage 
of this work. This work was partially financed by CNPq and FAPERGS, 
BRAZIL. 



References 

[1] Yu. L. Dokshitzer. Sov. Phys. JETP 46 (1977) 641; G. Altarelli and G. 
Parisi. Nucl. Phys. B126 (1977) 298; V. N. Gribov and L.N. Lipatov. 
Sov. J. Nucl. Phys 15 (1972) 438. 

[2] K. J. Eskola, V. J. Kolhinen, C. A. Salgado. Eur. Phys. J. C9 (1999) 
61; K.J. Eskola, V. J. Kolhinen, P. V. Ruuskanen. Nucl. Phys. B535 
(1998) 351. 



[3] A. H. Mueller, Nuc. Phys. B558 (1999) 285. 

[4] A.L. Ayala, V. P. Gongalves. Eur. Phys. J. C20 (2001) 343. 

[5] S. Liuti, R. Vogt, Phys. Rev. C51 (1995) 2244. 

[6] N. Hammon et al, Phys. Rev. C59 (1999) 2744. 

[7] M. B. Gay Ducati, V. P. Gongalves, Phys. Lett. B502 (2001) 92. 

[8] E. Iancu, A. Leonidov, L. McLerran, Nucl. Phys. A692 (2001) 583. 

[9] L. V. Gribov, E. M. Levin, M. G. Ryskin. Phys. Rep. 100 (1983) 1; A. 
H. Mueller, J. Qiu. Nucl. Phys. B268 (1986) 427; I. Balitsky, Nucl. 
Phys. B463 (1996) 99; A. L. Ayala, M. B. Gay Ducati and E. M. Levin. 
Nucl. Phys. B493 (1997) 305; Nucl. Phys. B511 (1998) 355; J. Jalilian- 
Marian, A. Kovner, H. Weigert, Phys. Rev. D59 (1999) 014015; Yu. 
Kovchegov, Phys. Rev. D60 (1999) 034008. 

[10] A. L. Ayala, M. B. Gay Ducati and E. M. Levin. Nucl. Phys. B493 
(1997) 305; Nucl. Phys. B511 (1998) 355. 

[11] S. Frixione et al., Nucl. Phys. B431 (1994) 453. 

[12] See, e. g., M. J. Leich et al, Nucl. Phys. A544 (1992) 197c. 

[13] R. Vogt, Phys. Rev. C61 (2000) 035203. 

[14] M. Gluck, E. Reya and A. Vogt. Z. Phys. C67 (1995) 433. 

[15] M. Gyulassy, L. McLerran, Phys. Rev. C56 (1997) 2219. 

[16] V. P. Gongalves, Phys. Lett. B495 (2000) 303. 

[17] L. McLerran, R. Venugopalan, Phys. Rev. D59 (1999) 094002. 

[18] F. Gelis and A. Peshier, Nucl. Phys. A697 (2002) 879. 

[19] A. Dumitru, L. McLerran, Nucl. Phys. A700 (2002) 492. 

[20] Conceptual Design for the Relativistic Heavy Ion Collider, BNL-52195, 
May, 1989; G. Baur et al, CMS Note 2000/060. 

[21] X-N. Wang, M. Gyulassy, Phys. Rev. Lett. 68 (1992) 1480. 



10 



Figure Captions 



Fig. |l|: Comparison between the ratio Rq with and without the high density 
effects for A=40 (Ca) and A=208 (Pb). The RHIC and LHC kinematic 
regions are also shown. 



Fig. |^: The exponent a [Eq. @] as a function of energy yfs for two 
nuclei (Ca and Pb) in pA processes. 



Fig. ||: The exponent a(xp) [Eq. (H)] as a function of Xp for p — Pb 
processes at RHIC (y/s = 350 GeV) and LHC (y/s = 8.8 TeV) (A=208) 
[Note the different inferior limit of the two graphs] . 



Fig. f|: The exponent (5 [Eq. (||)] as a function of energy yfs for Pb — Pb 
processes. 
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Figure 2: 
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